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STABILITY  CRITERIA  FOR  DIFFERENCE  APPROXIMATIONS  TO 
HYPERBOLIC  SYSTEMS,  AND  MULTIPLICATIVE  OF 
MATRIX  AND  OPERATOR  NORMS 


Principal  Investigator:  Moshe  Goldberg 


ABSTRACT 


Research  completed  under  Grant  AFOSR  -88-0175  by  Moshe 
Goldberg  consists  of  the  following  two  topics: 

(a)  Convenient  stability  criteria  for  difference  approximations  to 
hyperbolic  initial-boundary  value  problems. 

(b)  Multiplicativity  and  stability  of  matrix  and  operator  norms. 
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STABILITY  CRITERIA  FOR  DIFFERENCE  APPROXIMATIONS  TO 
HYPERBOLIC  SYSTEMS,  AND  MULTIPLICATIVE  OF 
MATRIX  AND  OPERATOR  NORMS 

Principal  Investigator:  Moshe  Goldberg 


1 .  Convenient  Stability  Criteria  for  Difference  Approximations  to 
Hyperbolic  Initial-Boundary  Value  Problems 


Consider  the  first  order  system  of  hyperbolic  partial  differential  equations 

8u(x,t)/dt  =  A3u(x,t)/5x  +  Bu(x,t)  +  f(x,t),  x>0,  t>0,  (1.1a) 

where  u(x,t)  =  (u(1)(x,t),  ....  u(n)(x,t))'  is  the  unknown  vector  (prime  denoting  the 

transpose),  f(x,t)  =  (f(1)(x,t) . f(n)(x,t))'  is  a  given  n-vector,  and  A  and  B  are 

fixed  n  x  n  matrices  such  that  A  is  diagonal  of  the  form 

a'  >  0,  a"  <  0,  (1.2) 


with  A1  and  A11  of  orders  kxk  and  (n  -  k)  x  (n  -  k),  respectively. 

The  solution  of  (1.1a)  is  uniquely  determined  if  we  prescribe  initial  values 

u(x,0)  =  u(x),  x  >  0,  (1.1b) 


and  boundary  conditions 


u"(0,t)  =  Su'(O.t)  +  g(t),  t  >  0, 


(lie) 
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where  S  is  a  fixed  (n  •  k)  x  k  coupling  matrix,  g(t)  a  given  (n  -  k)-vector,  and 


u'.(u"’ . A.  u"  =  <Ulk*” . U<"V 


(1.3) 


a  partition  of  u  into  its  outflow  and  inflow  components,  respectively, 
corresponding  to  the  partition  of  A  in  (1 .2). 

In  [GT7,  8]  E.  Tadmor  and  I  extended  our  results  in  [GT4-6]  and  obtained 
versatile,  easily  checkable  stability  criteria  for  a  wide  class  of  finite  difference 
approximations  to  the  above  initial-boundary  value  problem. 

More  specifically,  introducing  a  mesh  size  Ax  >  0,  At  >  0,  such  that  X  s  At/Ax 
is  a  constant,  and  using  the  notation  v(t)  =  v(jAx,t),  we  approximate  (1.1a)  by  a 
general,  basic  difference  scheme  -  explicit  or  implicit,  dissipative  or  not,  two- 
level  or  multilevel  -  of  the  form 


Q  , v (t  +  At)  =  T  Q  v (t  -  oAt)  +  Atb  (t),  j 

I  O  1  j 

o-0 


=  r,  r  +  1 . 


Q 

o 


i 

i  *  *  i 


A  EX, 


Ev 

i 


v  a  =  -1, ...,  s, 

i  + 1 


(1.4) 


where  the  n  x  n  coefficient  matrices  A  are  polynomials  in  XA  and  AtB,  and  the 
n-vectors  b  (t)  depend  on  f(x,t)  and  its  derivatives. 


The  difference  equations  in  (1.4)  have  a  unique  solution  v.(t  +  At)  if  we 
provide  initial  valuta 

v(oAt)  =  v(oAt),  o=0 . s,  j  =  0,1,2 .  (1.5) 

and  specify,  at  each  time  level  t  =  oAt,  o  =  s,  s  +  1 .  boundary  values 

v.(t  +  At),  j  =  0 . r  -  1.  Such  boundary  values  are  determined  by  boundary 

conditions  of  the  form 
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T^v. (t  +  At)  =  ^  T^v  (t  -  oAt)  +  Atd  (t),  j  =  0 . r  - 1 , 

J  a- 0 

m 

T(j)  V  r(j)  nT  „  i  „ 

T  =  /  .  C  E,  a  =  -1, ....  q, 

a  ■“  ot 

t-0 


(1.6a) 


where  the  n  x  n  matrices  C(J)  depend  on  A,  AtB  and  S,  and  the  n-vectors  d  (t) 

at  I 

are  functions  of  f(x,t),  g(t)  and  their  derivatives. 


Our  intention  was  to  interpret  the  difficult  and  often  stubborn  Gustafsson- 
Kreiss-Sundstrom  stability  criterion  in  [GKS]  in  order  to  obtain  convenient, 
simple  stability  criteria  for  approximation  (1.4)  -  (1.6a).  While  we  were  unable  to 
meet  this  goal  for  general  boundary  conditions  of  type  (1.6a),  we  managed  to 
achieve  rather  satisfactory  results  under  the  further  assumption  that,  in 

accordance  with  the  partition  of  A  in  (1 .2),  the  matrices  can  be  written  as 


where  for  B  =  0: 


ot 


_l(i) 

i  H  (i) 

c 

c 

ot 

ot 

C»(i) 

OT 

ot 

0  -6b) 


the  C1  (J>  are  independent  of  j, 

O  T 

(1.6c) 

the  C* (l)  are  diagonal, 

ot 

(1.6d) 

-inti)  . 

the  C  vanish  , 

o  t 

(1-6e) 

the  CM  (j)  vanish  for  t>0  orc>-1. 

o  t 

(16f) 

The  essence  of  (1.6c)-(1.6e)  is  that  for  B  =  0,  the  outflow  boundary 
conditions  are  translatory  (i.e.,  determined  at  all  boundary  points  by  the  same 
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coefficients),  separable  (i.e. ,  split  into  independent  scalar  conditions  for  the 
different  outflow  unknowns),  and  independent  of  outflow  values. 
Assumption  (l.6f)  implies  that  for  B  =  0,  the  inflow  values  at  the  boundary 
depend  essentially  on  the  outflow. 

It  should  be  pointed  out  that  our  outflow  boundary  conditions  are  quite 
general,  despite  the  apparent  restrictions  in  (1.6c)-(1.6e).  Indeed,  (1.6c)  is  not 
much  of  a  restriction,  since  in  practice  the  outflow  boundary  conditions  are  often 
translatory.  In  particular,  if  the  numerical  boundary  consists  of  a  single  point, 
then  the  boundary  conditions  are  translatory  by  definition,  so  (1.6c)  holds 
automatically.  The  restrictions  in  (1.6d)  and  (l.6e)  pose  no  great  difficulties 
either,  since  they  are  satisfied  by  all  reasonable  boundary  conditions,  where  for 

B  _  0  the  usually  reduce  to  polynomials  in  the  diagonal  block  A1,  and  the 

rm  (i) 

uot  vanish. 

We  realize  that  in  view  of  the  restriction  in  (1.6f)  our  inflow  boundary 
conditions  are  not  quite  as  general  as  the  outflow  ones.  They  can,  however,  be 
constructed  to  any  degree  of  accuracy  (see  [GT2});  and  they  can  be 
conveniently  extended  as  shown  in  [GT8].  We  note  that  if  the  boundary  consists 
of  a  single  point,  then  such  conditions  can  be  achieved  in  a  trivial  manner, 
simply  by  duplicating  the  analytic  condition  (1.1c),  i.e., 

v|J(t-KAt)  =  SVg(t  +  At)  +  g(t  +  At). 


Throughout  our  work  we  assume,  of  course,  that  the  basic  scheme  (1.4)  is 
stable  for  the  pure  Cauchy  problem,  and  that  the  additional  assumptions  which 
guarantee  the  validity  of  the  Gustafsson-Kreiss-Sundstrom  theory  in  [GKS], 
hold. 

The  first  step  in  our  analysis  was  to  reduce  the  above  stability  question  to 
that  of  a  scalar,  homogeneous  problem.  This  is  obtained  by  considering  the 
outflow  scalar  equation 

9u(x,  t)/8t  =  adu(x,  t)/3x,  x>0,  t  >  0,  a  =  constant  >  0,  (1.7) 


7 


for  which  the  basic  scheme  (1.4)  reduces  to  the  homogeneous  scheme 


S 

Q  .v  (t  +  At)  =  A  Qv(t-oAt) 

I  O  ] 

o»0 


Q 

a 


t  *  -  r 


a  =  -1 . s, 


(1.8a) 


and  the  boundary  conditions  (1 .6)  reduce  to  translatory  conditions  of  the  form 


T (t  +  At) 


i 


T  v  (t  -  aAt) 

o  i 


T 

<3 


-I 


c 

O  T 


O  =  -1 . q, 


(1.8b) 


where  a  and  c  are  scalar  coefficients. 

<J  T  Ot 

Referring  to  (1 .8)  as  the  basic  approximation,  we  proved: 

Theorem  1.1  [GT7,  8],  Approximation  (1.4)-(1.6)  is  stable  if  and  only  if 
the  basic  approximation  (1.8)  is  stable  for  every  eigenvalue  a  >  0  of  A1. 
That  is,  approximation  (1.4)-(1.6)  is  stable  if  and  only  if  the  scalar  outflow 
components  of  its  principal  part  are  all  stable. 

This  reduction  theorem  implies  that  from  now  on  we  may  restrict  our  stability 
study  to  the  basic  approximation  (1.8). 

In  order  to  introduce  our  stability  criteria  for  the  basic  approximation,  we  use 
the  coefficients  of  the  basic  scheme  (1.8a)  to  define  the  basic  characteristic 
function 


P(Z,K) 


t  » -  r 


s 


a  ■  0 
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Similarly,  using  the  coefficients  of  the  boundary  conditions  in  (1 .8b)  we  define 
the  boundary  characteristic  function 


R(z,k) 


t 


-a- 1 

C  Z 

ot 


T 

K  . 


Now  putting 


Q{z,k)  =  j  P(z,ic)  |  +  j  R(z,k)  I, 


we  proved: 

Theorem  1.2  [GT7,  8].  The  basic  approximation  (1 .8)  is  stable  if  either 


3P(z,k)  5P(z,k)  I  <  Q 

5z  5k  '  z=k«-i 

or 

Q(z  =  -1 ,  k  =  -1)  >  0  ; 

and  in  addition 


(1.10a) 


(1.10b) 


Q(z,k)  >0  for  all  |z|  —  |k|  =  1 ,  k  *  1,  (z,k)  *  (-1  ,-1),  (1 .10c) 

Q{z,  K=1)  >  0  for  all  |z|  =1,  z  *  1,  (1.1  Od) 

Q(z,k)  >  0  for  all  |z|  >1,0  <  |  k|  <  1.  (1.1  Oe) 

The  advantage  of  this  setting  of  Theorem  1.2  is  clarified  by  the  following 
lemma,  in  which  we  provide  helpful  sufficient  conditions  for  each  of  the  four 
inequalities  in  (1.10b  -  e)  to  hold: 
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Lemma  1.3  [GT7,  8]. 

(i)  Inequalities  (I.IOb.c)  hold  if  either  the  basic  scheme  (1.8a)  or  the 
boundary  conditions  (1 .8b)  are  dissipative. 

(ii)  Inequality  (1.1  Od)  holds  if  any  of  the  following  is  satisfied: 

(a)  The  basic  scheme  is  two-level. 

(b)  The  basic  scheme  is  three-level  and 

Q(z  =  -1  ,k  =  1)  >  0.  (1.11) 

(c)  The  boundary  conditions  are  at  most  two-level  and  at  least 
zero-order  accurate  as  an  approximation  to  equation  (1.7). 

(d)  The  boundary  conditions  are  three-level,  at  least  zero-order 
accurate  as  an  approximation  to  ( 1.7),  and  (1.11)  /s  satisfied. 

(iii)  Inequality  (l.lOe)  holds  if  the  boundary  conditions  fulfill  the  von 
Neumann  condition,  and  are  either  explicit  or  satisfy 


T,(k) 


*  0 


V  0<|K|<1. 


The  stability  criteria  obtained  in  Theorem  1.2  depend  both  on  the  basic 
scheme  and  the  boundary  conditions,  but  not  on  the  intricate  and  often 
complicated  interaction  between  the  two.  Consequently,  Theorem  1.2,  aided  by 
Lemma  1.3,  provide  in  many  cases  a  convenient  alternative  to  the  celebrated 
stability  criteria  of  Kreiss  [K2]  and  of  Gustafsson,  Kreiss  and  Sundstrom  [GKS]. 

Having  the  new  criteria,  we  easily  established  stability  for  a  host  of  examples 
that  incorporate  and  generalize  most  of  the  cases  studied  in  recent  literature; 
e.g.,  [Gol ,  Go2,  GGT,  GKS,  GO,  GTi ,  GT2,  GT4-6,  K1 ,  KO,  01,  Osh,  Ski ,  Sk2, 
SK,  Ta,  Th.Trl],  To  mention  just  a  few  of  our  examples,  we  obtained  stability  for: 

(a)  Any  stable  basic  scheme,  with  boundary  conditions  generated  by 
either  the  explici4  or  implicit  one-sided  Euler  schemes. 

(b)  Any  stable  two-level  basic  scheme,  with  boundary  conditions 
generated  by  either  horizontal  extrapolation  or  by  the  one-sided 
three-level  Euler  scheme. 
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(c)  Any  stable  dissipative  basic  scheme,  with  boundary  conditions 
generated  by  oblique  extrapolation  or  by  the  Box  scheme. 

(d)  The  Crank-Nicolson,  Backward-Euler,  Leap-Frog  and  Lax-Friedrichs 
schemes  (all  nondissipative),  with  boundary  conditions  generated  by 
either  oblique  extrapolation  or  by  the  one-sided  Weighted  Euler 
scheme. 


We  drew  great  satisfaction  from  the  fact  that  our  theory  and  examples  in 
[GT4-6]  were  used  already  by  a  number  of  authors,  including  Berger  [Bm], 
LeVeque  [Le],  South,  Hafez  and  Gottlieb  [SHG],  Thune  fThl ,  2],  Trefethen  [Trl , 
2],  and  Yee  [Y],  Thune,  in  his  effort  to  provide  a  software  package  for  stability 
analysis  of  finite  difference  approximations  to  hyperbolic  initial-boundary  value 
problems,  says  in  [Th  1  ]:  "...Another  approach  has  been  to  derive  new  criteria, 
based  on  the  Gustafsson-Kreiss-Sundstrom  theory  but  more  convenient  for 
practical  use...  The  most  far-reaching  work  along  these  lines  has  been  made  by 
Goldberg  and  Tadmor  [GT1,  2,  4]  ..."  In  [Th2]  he  says:  "A  conceptually  different 
approach  has  been  suggested  by  Goldberg  and  Tadmor  [GT6].  They  have 
treated  a  fairly  general  class  of  difference  approximations  of  problems  in  one 
space  dimension.  By  exploring  the  properties  of  this  class  they  have  been  able 
to  ...  simplify  the  stability  analysis  considerably".  Talking  about  the  ultimate 
goal  of  his  project,  Thune  [Th2]  adds:  "...the  black  box  design  ought  to  be 
abandoned.  A  flexible  tool  box  design  would  be  preferable.  For  example,  the 
environment  shoud  include  tools  for  checking  the  convenient  stability  conditions 
of  Goldberg  and  Tadmor". 

We  were  also  pleased  to  learn  that  part  of  our  theory  in  [GT7]  was  taught 
already  in  several  institutions  including  UCLA,  NYU,  and  the  University  of 
Paris  VI. 


2.  Norms,  Seminorms  and  Multiplicativity  Factors 

Let  A  be  an  algebra  over  the  complex  field  C.  As  usual,  a  real-valued 
function 


it 


S  :  A  -»R 


is  called  a  seminorm  if  for  all  x,  y  e  a  and  X  e  C: 

S(x)  >  0, 

S(Xx)  =  |  X  |  •  S(x), 

S(x  +  y)  <  S(x)  +  S(y). 

If  in  addition,  S  is  positive  definite,  i.e., 

S(x)  >0,  x  *  0, 

then  S  is  a  norm.  We  call  a  seminorm  S  proper  if  S  *  o  and  S(x)  =  0  for  some 
x  *  0.  Finally,  we  say  that  S  is  submultiplicative  (or  simply,  multiplicative  )  if 

S(xy)  <  S(x)S(y),  Vx.ye  A. 

Of  special  interest  are,  of  course,  operator  algebras.  Here,  we  have  a 
normed  vector  space  V  over  C  and  an  algebra  «(V)  of  bounded  linear 
operators  on  V. 

The  first  example  that  comes  to  mind  of  a  multiplicative  norm  on  an  operator 
algebra  is  the  ordinary  operator  norm. 


II  A  ||  =  sup  {  I  Ax  |  :  x  e  V,  |  x  |  =  1  }, 


(2.1) 


where  |  |  is  the  vector  norm  on  V. 

If  V  is  a  Hilbert  space,  then  a  well  known  example  of  a  nonmultiplicative 
norm  on  «(V)  is  the  numerical  radius  (e.g.,  [Be,  Go3,  GT3,  Ha,  P] ): 

r(A)  =  sup  {  |  (Ax,  x)  |  :  x  e  V,  |  x  |  =  (x,  x)1/2  =  1  }  ,  (2.2) 

which  plays  an  important  role  in  stabil,  /  analysis  of  finite  difference  schemes  for 
multi-space-dimensional  hyperbolic  initial-value  problems  [GT3,  Li,  LW,  Tu], 
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Another  example  of  considerable  interest  is  the  l  norm,  1  <  p  <  «,  defined 
on  Cn  x  ,  the  algebra  ofnxn  complex  matrices: 


A  =  (a.-)  e  C 

•I  n  x  n 


(2.3) 


Ostrowski  [Ost]  has  shown  that  this  norm  is  multiplicative  if  and  only  if  1  <  p  <  2. 

Given  a  seminorm  S  on  an  arbitrary  algebra  A,  and  a  fixed  constant  p  >  0, 
then  obviously  =nS  is  a  seminorm  too.  Clearly,  may  or  may  not  be 
multiplicative.  If  it  is,  we  call  p  a  mjltiplicativity  factor  for  S.  That  is,  p  is  a 
multiplicativity  factor  for  S  if  and  only  if 

S(xy)  <  pS(x)S(y),  Vx,y  e  A. 

Evidently,  if  p0  is  a  multiplicativity  factor  of  S,  then  so  is  any  p  with  p  >  p0. 
Thus,  having  a  seminorm  S,  the  question  is  whether  S  has  multiplicativity 
factors;  and  if  so,  is  there  a  best  (least)  one? 

This  question  can  be  easily  answered  as  follows: 


Theorem  2.1  ([GS1],  [AG1  ]).  Let  A  be  an  algebra,  and  let  S  *  0  be  a 
seminorm  on  A.  Then: 

(a)  S  has  multiplicativity  factors  if  and  only  if  Ker  S  is  an  ideal  in  a  and 

P  inf  =  sup{S(xy )  :  x,y  e  A,  S(x)=S(y)  =  1}<  «  .  (2.4) 

(b)  If  S  has  multiplicativitiy  factors  and  Pjnf  >0,  then  p  is  a 

multiplicativity  factor  if  and  only  if  p  >  pjn}. 

(c)  If  S  has  multiplicativitiy  factors  and  Pinf  =  o,  then  p  is  a 

multiplicativity  factor  if  and  only  if  p  >  0- 
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If  a  is  finite  dimensional  and  S  is  a  norm,  then  Ker  S  =  {0}  is  an  ideal  in  a, 

and  a  standard  compactness  argument  implies  that  hnf  in  (2.4)  is  finite;  so  by 
Theorem  2.1 ,  S  has  multiplicativity  factors.  This  is  not  always  the  case  if  S  is  a 
proper  seminorm. 

In  the  infinite  dimensional  case,  it  was  shown  by  Straus,  Arens  and  myself 
[GS2,  AG1],  that  norms  as  well  as  proper  seminorms  may  or  may  not  have 
mutiplicativity  factors. 

While  Theorem  1.1  seems  to  settle  the  question  of  characterizing 
multiplicativity  factors,  the  quantity  pjnf  in  (2.4)  is  often  difficult  to  compute.  A 

more  practical  approach  toward  checking  whether  a  given  constant  p  >  0  is  the 
best  (least)  multiplicativity  factor  for  a  given  seminorm  S,  is  by  verifying  that 

S(xy  )  <  pS(  x  )S(  y  )  Vx.ye 

with  equality  for  some  x0  and  y0  for  which  S(x0)  *  0,  S(y0)  *  0. 

Using  this  observation,  Holbrook  [Ho],  and  independently  Straus  and  myself 
[GS1],  showed  that  if  V  is  a  Hilbert  space  over  C  of  dimension  at  least  2,  and  r  is 

the  numerical  radius  in  (2.2),  then  the  best  multiplicativity  factor  for  r  is  p  =4. 

Similarly,  Maitre  [M],  and  Straus  and  I  [GS3],  showed  that  the  best 
multiplicativity  factor  for  the  /p  norm  on  Cn  x  n  defined  in  (2.3),  is: 

j  1  ,  1<p<2 

\  n1  2/p  ,  2  <  p  <  oo  . 

Often,  when  the  least  multiplicativity  factor  remains  unknown,  one  may  try  to 
obtain  multiplicativity  factors  through  the  following  version  of  a  theorem  by 
Gastinel. 

Theorem  2.2  [Ga,  GS2,  AG1].  Let  S,  T  be  seminorms  on  an  algebra 
A,  Let  T  be  multiplicative,  and  let  z  >  o  >  0  be  constants  such  that 

oT(x  )  <S(x)  £tT(x),  Vxe  a. 

Then  any  p  >  t  /  a  2  in  a  multiplicative  factor  for  S. 
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This  result  was  utilized  by  Straus  and  myself  [GS2,  4]  to  obtain 
multiplicativity  factors  for  certain  generalizations  of  the  numerical  radius,  called 
C-numerical  radii. 

In  [AG1]  Arens  and  I  investigated  multiplicativity  factors  in  terms  of  the 
kernels  of  S.  In  particular  we  proved: 

Theorem  2.3  [AG1J.  If  a  is  a  simple  algebra  then  there  are  no 
multiplicative  proper  seminorms  on  a. 


Since  Cnxn  is  simple  (e.g.  [BM]  ),  we  immediately  obtain  from  Theorem  2.3 
the  following  result  that  was  directly  proved  by  Straus  and  myself  in  [GS1]: 

Theorem  2.4  [AG1].  There  are  no  multiplicative  proper  seminorms  on 

C 

n  x  n 

Arens  and  I  [AG1,  2]  specialized  our  study  to  function  algebras  and  to 
seminorms  generated  by  the  sup  norm,  where  we  gave  the  following  three 
characterizations  of  multiplicativity  factors: 

Theorem  2.5  [AG1,  2].  Let  T  be  a  set  and  let  a  be  the  algebra  of 
bounded  functions 


x:T  ->  C, 

with  the  usual  multiplication 

xy(  t )  =  x(  t  )y(  t );  x,  y  e  a:  t  e  T. 

For  a  fixed  element  c,  0  *  c  e  a,  define  the  seminorm 

S  (x)  *  sup  |c(t)x(t)|  .  (2.5) 

te  T 


15 


Then: 


(a)  Sc  has  muttipUcativity  factors  if  and  only  if 

e  3 inf{|c(  t  )j  :  t  e  T,  c(t  )*0}  >0. 

(b)  Ife  >  0,  then  p  >  0  is  a  muttipUcativity  factor  for  Sc  if  and  only  if 
p  >  e 

Theorem  2.6  [AG1,  2].  Let  T  be  a  topological  space  and  let  A  be  the 
algebra  of  bounded  continuous  functions 

x:  T  ->  C  . 

For  a  fixed  c,  0  *  c  e  a,  let  Sc  be  the  seminorm  in  (2.5).  Then  conclusions 
(a)  and(b)  of  Theorem  2.5  hold. 

Theorem  2.7  [AG1,  2).  Let  T,  a,  c,  and  Sc  be  as  in  Theorem  2.6,  and 
suppose  T  is  connected.  Then 

(a)  The  following  are  equivalent: 

(i)  Sc  has  the  muttipUcativity  factors. 

(ii)  5einf{|c(t  )|:  t  €  T}  >  0. 

(iii)  c  is  invertible. 

(iv)  Sc  is  a  norm  on  a. 

(b)  If  5  >  0  then  p  >  0  is  a  multiplicativity  factor  for  Sc  if  and  only  if 
4  >  5'1  . 

For  example,  consider  i°°,  the  algebra  of  bounded  sequences  x  =  {  t 
over  C,  with  the  usual  Hadamard  multiplication 


xy  =  fe},  x  =  {U.y  =  {TUe  r. 

11  j  J 
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Fix  an  element  c  =  {y}  e  l°°  c  *  0,  and  define  the  seminorm, 

i  9 


sc(x)  =  sup|  Yj^jl.  X  =  (4j}  €  r . 
Obviously,  Sc  is  a  norm  on  i°°  if  and  only  if 


Yj  *  0,  j. 1.2.3 . 

Otherwise  Sc  is  a  proper  seminorm. 

By  Theorem  2.5  (here  T  =  2+  =  {1,2,3,  ... }),  Sc  has  multiplicativity  factors  if 
and  only  if 


e  *  inf  I  Y  j  1  >  0  , 

Y  ,  *0 


and  if  e  >  0  then  the  best  (least)  multiplicativity  factor  for  Sc  is  pmin  =  e'1  . 

The  four  simple  selections 

V  1.  j  =  1 .2,3,  ... 

yj  =  j'\  j  =  i ,2,3, ... 
y,  =0;  7j  *1.  j  =2,3,4,... 

Y:  =  Y j  =  i'\  i  =2,3,4,  ... 

show,  as  indicated  before,  that  in  the  infinite  dimensional  case,  both  norms  and 
proper  seminorms  may  or  may  not  have  multiplicativity  factors. 
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Next,  let  S  be  a  seminorm  on  Cn  and  consider  the  real-valued  function  S+ 
defined  by 


S+(x)  =  S(x+),  xeC", 

where 

X+  =  (1^1)  €  R" 

is  the  vector  obtained  by  taking  the  absolute  values  of  the  components  of 

As  in  [Go9],  we  call  S  quasimonotonic  if 

x,  y  €  Rn  with  0  <  x  <  y  implies  S(x)  <  S(y), 
where  the  inequalities  0  <  x  <  y  are  construed  componentwise. 

With  this  definition  we  can  prove: 

Theorem  2.8  [Go9].  Let  S  be  a  seminorm  (norm)  on  Cn.  Then  S+  is  a 
seminorm  (norm)  on  Cn  if  and  only  if  S  is  quasimonotonic. 

For  example,  it  was  shown  in  [Go8]  that  the  numerical  radius  r  in  (2.2)  is 
quasimonotonic  on  Cn  x  n.  Since  r  is  a  norm,  so  is  r+  by  Theorem  2.8. 

In  order  to  discuss  multiplicativity  factors  for  S+  let  us  assume  that  Cn  has 
been  given  a  structure  of  an  algebra  over  C.  This  can  be  done,  for  instance,  by 
taking  Cn  with  the  usual  Hadamard  multiplication 

xy  -  (SjTij),  X  =  ({;.),  y  =  (t|.  )  e  C", 

or  Cnx  n  with  usual  matrix  multiplication. 

We  can  now  prove  the  following  result  that  associates  multiplicativity  factors 
of  S  with  those  of  S+: 
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Theorem  2.9  [Go9].  Let  S  be  a  quasimonotonic  seminorm  on  Cn  with  a 
muttiplicativity  factor  p,  and  suppose  that  Cn  has  been  given  a  structure  of  an 
algebra  such  that 

(xy)+  <  x+y+  V  x,  y  e  C"  .  (2.6) 


Then: 

(a)  p  is  a  muttiplicativity  factor  for  S+. 

(b)  If  S  has  a  least  muttiplicativity  factor,  then  so  does  S+,  and  these 
least  factors  satisfy 


(S)  2  Hmin  (S*).  (2.7) 

Let  us  point  out  that  condition  (2.6)  holds  tor  all  common  multiplication  rules 
on  Cn,  including  Hadamard's  multiplication  on  Cn  and  the  standard  matrix 
multiplication  on  Cnxn  where  we  have 

(xy)  =  x  y  Vx.yeC 

and 

(AB)+  <  A+B+  VA,  BeC  , 

n  x  n 


respectively. 

We  further  remark  that  the  numerical  radius  in  (2.2)  satisfies 

Hnm  l'>  -  =  4' 

so  equality  in  (2.7)  is  possible.  It  was  shown,  however,  in  [Go8]  that  in  general 
the  ratio 


u  (S)/u  (S+) 

^mm  '  '  ^mtn  '  ' 


can  be  arbitrarily  large. 
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The  above  concepts  of  multiplicative  and  multiplied, factors  can  be 
extended  as  follows. 

De.ini.ion,  Le,  U  VttndWbe 

:;:r:  *■  —  -<•  ^  r:  * 

Z  sem, norms  on  s,  and  *  respectively,  and  p  >  0  is  a  constant  such  tha, 
S^AB)  <  hS2(A)S3{B)  VAe%  Be  <83, 
then  we  say  tha,  p  is  a  mu»iPlicaMy  factor  for  S,  with  respect  to  S2  and  S, 

jaassass  ausss*' 

r(AB)  <  2r(A)j|  B  !|  V  A,  B  e  ®(V)t 

with  equality  for  certain  operators  A  *  0  and  B  *  0.  Thus,  n-  2  is  the  best 
multiplicativity  factor  lor  r  with  respect  to  r  and  II  II- 

demo^r^^ 

1  <  p  <  the  1  norm  of  an  m  x  n  matrix  A  «  (a,)  -  „«■ 


I  m 

-  IXi«^ 

1  *  4 


i « 1  j  *  1 


(2.8) 


Defining 


P-P  - 


We  have: 


Theorem  2.10  [Go4],  Let  p,  q,  r  satisfy  1  <  p,  q,  r  £  and  let  q’  be  the 
conjugate  of  q  (/'.©.,  1/q  +  1/q'  =  1).  Then  the  best  multiplicativity  factor  for 
the  l  norm  on  C  with  respect  to  the  l  norm  on  C  .  and  the  l  norm 

p  mxn  q  mxk  r 

°nC**„is 


(1  =  XM(m)  y  n)  Xq.,(k). 

That  is,  for  all  A  e  C  .  and  Be  C, 

mxk  k  x  n 

I  AB  |p  <  XM(m)  yn)  Xq,(k)  |  A  |q  |  B  |r  ,  (2.9) 

where  this  inequality  is  sharp. 

Theorem  2.10  has  quite  a  few  applications.  For  example  (see  [Go6,  7]), 
taking  (2.9)  with  m  =  n  =  1,  we  get  an  upper  bound  for  the  standard  inner 
product  (x,  y)  on  Cn  in  terms  of  |  x  |  and  |  y  |p;  and  if  we  further  set  r  =  q’  we 
obtain  the  classical  Holder  inequality. 

Another  application  of  Theorem  2.10  concerns  the  ordinary  !p  operator-norm 


A  ||p  =  max  {  |  Ax  |p  :  x  e  C".  |  x  |p  =  1  } , 


(2.10) 


for  which  we  obtain: 


Theorem  2.11  [Go7],  Let  p,  q,  r  satisfy  1  <  p,  q.  r,  <  ».  Then  for  all 
A  €  C  . ,  B  e  C.  , 

mxk  kxn 

II  AB  ||p  <  Xw(m)  xqp(k)  yk)  y  n)  ||  A  ||q  ||  B  ||f  , 


where  the  inequality  is  sharp  if  either  q<p<r  or  r  <  p  <  q. 

A  third  consequence  of  (2.9)  describes  the  equivalence  relations  between 
the  norms  in  (2.8)  and  (2.9): 
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Theorem  2.12  [Go5].  Let  p,  q  satisfy  1  ^  p,  q  <  «,  and  let  q'  be  the 
conjugate  of  q.  Then  for  all  A  e  Cm  x  n, 

lAlp^^pq(mn)!Alq- 

HAHp^Wm)Vn)l|A|lq’ 

ilAHp^Wm)Vn)!A|q- 

|A)p<{mn)1/p||A)|q, 

where  the  first  three  inequalities  are  sharp. 
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22.  Invited  speaker,  The  1984  Haifa  Conference  on  Matrix  Theory,  Technion  - 
Israel  Institute  of  Technology  and  the  University  of  Haifa,  Haifa,  Israel, 
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31.  Invited  speaker  and  Session  Chairman,  First  International  Conference  on 
Industrial  and  Applied  Mathematics,  Minisymposium  on  Finite  Difference 
Approximations  to  Hyperbolic  Initial-boundary  Value  Problems,  Paris, 
June-July  1987,  title:  "Convenient  stability  criteria  for  difference 
approximations  of  hyperbolic  initial-boundary  value  problems". 

32.  Invited  speaker,  Meeting  on  Numerical  Problems  for  initial  and  Initial¬ 
boundary  Value  Problems,  Mathematics  Research  Institute,  Oberwolfach, 
West  Germany,  August  1987,  title:  "Stability  criteria  for  finite  difference 
approximations  to  hyperbolic  initial-boundary  value  problems". 

33.  Invited  speaker,  The  Fourth  Haifa  Matrix  Theory  Conference,  Technion  - 
Israel  Institute  of  Technology,  Haifa,  Israel,  January  1988,  title:  "On 
monotone  and  semi-monotone  matrix  functions". 

34.  Speaker  and  Session  Chairman,  Second  International  Conference  on 
Hyperbolic  Problems,  RWTH  Aachen,  Aachen,  West  Germany,  March 

1988,  title:  "Co p'»  tent  stability  criteria  for  difference  approximations  of 
hyperbolic  iniha1  1  uundary  value  problems". 

35.  Invited  speaker  and  Session  Chairman,  The  1988  Annual  Meeting  of  the 
Israel  Mathematical  Union,  Tel  Aviv  University,  Tel  Aviv,  Israel,  March 
1958,  title:  "Simple  stability  criteria  for  difference  approximations  of 
hyperbolic  initiai  joundary  value  problems". 

36.  Invited  speaker  and  Session  Chairman,  The  Fifth  Haifa  Matrix  Theory 
Conference,  Technion-lsrael  Institute  of  Technology,  Haifa,  Israel,  January 

1989,  title:  "Multiplicativity  factors  for  seminorms". 

37.  Invited  speaker,  Technion  Linear  Algebra  Summer  Workshop,  Technion  - 
Israel  Institute  of  Technology,  Haifa,  Israel,  June  1989,  title:  "On  a  class  of 
absolute  norms". 

38.  Invited  speaker  and  Session  Chairman,  Inaugural  Conference  of  the 
International  Linear  Algebra  Society,  Brigham  Young  University,  Provo, 
Utah,  August  1989,  titie:  "Norms,  seminorms  and  multiplicativity  factors". 
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2.  A  note  comparing  the  root  condition  and  the  resolvent  condition  (with  W.L. 
Miranker),  Information  Sciences  4  (1972),  285-288. 

3.  A  note  on  the  stability  of  an  iterative  finite-difference  method  for  hyperbolic 
systems,  Mathematics  of  Computation  27  (1973),  41-44. 

4.  Stable  approximations  for  hyperbolic  systems  with  moving  internal 
boundary  conditions  (with  S.  Abarbanel),  Mathematics  of  Computation  28 
(1974),  413-447. 

5.  Stable  schemes  for  hyperbolic  systems  with  moving  internal  boundaries 
(with  S.  Abarbanel),  in  "Computational  Mechanics",  edited  by  J.T.  Oden, 
Texas  Institute  for  Computational  Mechanics  (TICOM),  1974,  469-478. 

6.  On  matrices  having  equal  spectral  radius  and  spectral  norm  (with 
G.  Zwas),  Linear  Algebra  and  Its  Applications  8  (1974),  427-434. 

7.  The  numerical  radius  and  spectral  matrices  (with  E.  Tadmor  and  G.  Zwas), 
Linear  and  Multilinear  Algebra  2  (1975),  317-326. 
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8.  Numerical  radius  of  positive  matrices  (with  E.  Tadmor  and  G.  Zwas),  Linear 
Algebra  and  Its  Applications  12  (1975),  209-214. 

9.  On  inscribed  circumscribed  conics  (with  G.  Zwas),  Elemente  der 
Mathematik  31  (1976),  36-38. 

10.  Inclusion  relations  between  certain  sets  of  matrices  (with  G.  Zwas),  Linear 
and  Multilinear  Algebra  4  (1976),  55-60. 

11.  A  test  problem  for  numerical  schemes  for  nonlinear  hyperbolic  equations 
(with  S.  Abarbanel),  Computer  Methods  in  Applied  Mechanics  and 
Engineering  8  (1976),  331-334. 

12.  Inclusion  relations  involving  k-numerical  ranges  (with  E.G.  Straus),  Linear 
Algebra  and  Its  Applications  15  (1976),  261-270. 

13.  On  characterizations  and  integrals  of  generalized  numerical  ranges  (with 
E.G.  Straus),  Pacific  Journal  of  Mathematics  69  (1977),  45-54. 

14.  On  a  boundary  extrapolation  theorem  by  Kreiss,  Mathematics  of 
Computation  31  (1977),  469-477. 

15.  Elementary  inclusion  relations  for  generalized  numerical  ranges  (with  E.G. 
Straus),  Linear  Algebra  and  Its  Applications  18  (1977),  1-24. 

16.  On  a  theorem  by  Mirman  (with  E.G.  Straus),  Linear  and  Multilinear  Algebra 
5  (1977),  77-78. 

17.  On  boundary  extrapolation  and  dissipative  schemes  for  hyperbolic 
problems,  Proceedings  of  the  1977  U.S.  Army  Numerical  Analysis  and 
Computer  Conference,  ARO  Report  77-3  (1977),  157-164. 

18.  Scheme-independent  stability  criteria  for  difference  approximations  of 
hyperbolic  initial-boundary  value  problems.  I,  (with  E.  Tadmor), 
Mathematics  of  Computation  32  (1978),  1097-1107. 

19.  Norm  properties  of  C-numerical  radii  (with  E.G.  Straus),  Linear  Algebra 
and  Its  Applications  24  (1979),  1 13-132. 
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20.  On  certain  finite  dimensional  numerical  ranges  and  numerical  radii,  Linear 
and  Multilinear  Algebra  7  (1979),  329-342. 

21.  Combinatorial  inequalities,  matrix  norms,  and  generalized  numerical  radii 
(with  E.G.  Straus),  in  "General  Inequalities  2",  edited  by  E.F.  Beckenbach, 
Birkhauser  Verlag,  Basel,  1980,  37-46. 

22.  Scheme-independent  stability  criteria  for  difference  approximations  of 
hyperbolic  initial-boundary  value  problems.  II,  (with  E.  Tadmor), 
Mathematics  of  Comput  tion  36  (1981),  603-626. 

23.  On  the  numerical  <adius  and  its  applications  (with  E.  Tadmor),  Linear 
Algebra  and  Its  Applications  42  (1982),  263-284. 
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PROBLEMS  IN  APPLIED  AND  COMPUTATIONAL 
MATRIX  AND  OPERATOR  THEORY 


Principal  Investigator:  Marvin  Marcus 


ABSTRACT 


Research  completed  under  Grant  AFOSR-88-01 75  by  Marvin 
Marcus  during  the  period  5/1/88  -  11/30/90  consists  of  the 
following  topics: 

(a)  Hadamard  Products  and  Powers 

(b)  Inequalities  for  Tensors 

(c)  Inequalities  for  Generalized  Matrix  Functions 

(d)  Inequalities  for  Eigenvalues  and  Singular  Values 

(e)  Distance  Matrices 

(f)  Numerical  Range 

(g)  Determinants  of  Sums 

In  the  following  report  papers  are  listed  that  were  partially  or 
entirely  completed  during  the  reporting  period. 
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EIGENVALUE  AND  SINGULAR  VALUE  INEQUALITIES,  TENSOR 
SPACES,  AND  THE  NUMERICAL  RANGE 


Principal  Investigator:  Marvin  Marcus 


Each  of  the  sections  that  follow  begin  with  the  current  disposition  of  the  research  paper  named  in 
the  section  heading.  A  description  of  the  research  and  appropriate  references  follow  each 
section  heading. 

1  .  Hadamard  Square  Roots 

This  paper  is  currently  in  press  in  the  SIAM  Journal  of  Matrix  Analysis.  The  proof  sheets  have 
been  corrected  and  returned  to  the  editor  during  the  period  of  this  report. 

If  A  is  an  n-square  positive  semi-definite  hermitian  matrix  of  rank  1  then  the  Hadamard  square  root 
of  A  is  the  n-square  matrix  obtained  by  replacing  each  entry  of  A  by  the  principal  value  of  its  square 
root.  It  is  proved  that  if  A  has  no  zero  or  negative  entries  then  the  Hadamard  square  root  has  odd 
rank  and  all  odd  ranks  are  possible. 


REFERENCES 

[1]  P.  R.  Halmos,  Finite-Dimensional  Vector  Spaces.  2nd.  ed.,  D.  Van  Nostrand  Company, 
Inc.,  Princeton,  New  Jersey,  1958. 

[2]  R.  A.  Horn,  The  Hadamard  Product,  Proceedings  of  Symposia  in  Applied  Mathematics, 
American  Mathematical  Society,  to  appear. 

[3]  R.  A.  Horn  and  C.  R.  Johnson,  Matrix  Analysis,  Cambridge  University  Press,  New  York. 
1985. 

[4]  M.  Marcus  and  N.  Khan,  A  note  on  the  Hadamard  product,  Canad.  Math.  Bull.,  2  (1959) 
pp.  81-83. 
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2 .  Multilinear  Methods  in  Linear  Algebra 

This  paper  is  currently  in  press  in  the  journal  Linear  Algebra  and  its  Applications.  It  is  a  written 
version  of  a  one  hour  invited  address  at  the  first  meeting  of  the  International  Linear  Algebra 
Society  held  at  Brigham  Young  University,  Provo  Utah,  8/12/89-  8/16/89. 

Several  classical  and  new  results  are  presented  in  which  multilinear  algebra  has  proven  to  be  an 
effective  tool.  Some  of  the  topics  covered  are:  Weyl’s  inequalities:  the  Hadamard  product: 
mappings  on  tensor  spaces;  Gram  matrices  in  tensor  spaces:  strong  non-singularity  and  the  LDU 
theorem. 
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81  Wielandt,  Helmut,  Topics  in  the  Analytic  Theory  of  Matrices,  Department  of  Mathematics, 
University  of  Wisconsin,  Madison,  1967. 

3 .  Experiments  with  Entrywlse  Square  Roots 

This  paper  is  currently  scheduled  to  appear  in  ;he  next  issue  of  The  Math  Works  Newsletter.  The 
Math  Works  Inc.  distributes  the  linear  algebra  package,  MATLAB. 

In  this  paper  MATLAB™  was  used  to  compute  the  entrywise  square  root  of  an  n-square  complex 
matrix  A  and  to  investigate  some  of  the  usual  invariants,  e.g.,  rank  (sqrt(A)).  A  particularly 
important  class  of  matrices  are  the  n-square  positive  semi-definite  hermitian  A  of  rank  1  (such 
matrices  are  the  heart  of  the  matter  for  the  normal  spectral  decomposition  theorem).  It  is  easy  to 
generate  random  complex  A  of  this  kind: 

x  =  2  *  rand  (1,n)  -  ones  (1,n); 
y  =  2  *  rand(1,n)  -  ones  (1,n); 
u  =  x  +  i  *  y; 

A  =  u' »  u; 

Of  course,  if  A  has  positive  entries  there  is  no  difficulty  in  confirming  that  sqrt(A)  has  rank  1 .  If  A 
has  a  negative  entry  then  sqrt(A)  need  not  be  hermitian,  nor  of  rank  1 : 
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In  fact,  there  are  A  >  0  (positive  semi-definite  hermitian)  with  positive  entries  for  which  sqrt(A)  is 
indefinite  [2,  p.  462]: 


100 
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1 

9 

100 

0 

81 

4 

0 

100 

16 

1 

81 

16 

100 

The  eig(A)  command  produces 

1 .0205354435388008+02 
9.771678596998484e+01 
1.832428791889355e+02 
1.6986790487199730+01 

whereas  eig(sqrt(A))  yields 

1.0335731941645718+01 

9.0275667247986940+00 

-1.8733739480494150-01 

2.0824038728360540+01 

The  MATLAB™  code  exhibited  above  can  easily  be  embedded  in  a  script  to  compute  and 
tabulate  rank(sqrt(A))  for  a  large  number  of  complex  n-square  A  >  0,  rank  (A)  =  1 .  What  is 
surprising  is  that  for  all  n  tested,  and  significantly  large  samples  (250  random  A),  rank(sqrt(A))  was 
computed  as  an  odd  integer  in  every  case.  These  experimental  results  led  to  the  conjecture  and 
proof  of  the  following  result. 

If  A  is  an  n-square,  rank  i ,  positive  semi-definite  hermitian  matrix  with  no  negative 
entries  then  rank(sqrt(A))  is  always  an  odd  integer.  Moreover,  if  v  is  any  odd  integer, 

1  <  v  <  n,  there  exists  an  n-square  A  >  0,  rank(A)  =  1 ,  with  no  zero  or  negative 
entries,  for  which  rank(sqrt(A))  =  v. 
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4.  Bessel's  Inequality  in  Tensor  Space 

This  paper  appeared  during  the  period  of  the  current  report:  Linear  and  Algebra 

(1988)  Vol.  23,  pp.  233-249. 

Let  A  be  an  n-square  complex  matrix  and  define  to  be  the  n!-square  matrix  h^se  entries  are 

n 

IT-  ao(i)  i(i) 1  where  a  ancl  T  njn  lexicographically  over  Sn.  If  A  is  positive  definite  hermitian  and  x  is 
a  unit  n!-tuple  then 


(\X,X)  *  det(A)  +  |  X  X(o)  I2  c(A) 


where  c(A)  is  the  largest  of  the  numbers  |  a  /  a°  ,  j  =  1 . n,  and  the  summation  is  over  o 

e  Sn  For  n  =  3,  if  A  is  not  permutation  similar  to  a  direct  sum  and  ^  is  a  unit  n!-tuple  then 

(:Tax,x)  =  det(A)  iff  x  is  a  multiple  of  the  alternating  character  The  relationships  among  recent 
results  of  Bapat  and  Sunder,  Chollet,  and  Gregorac  and  Hentzel  are  also  discussed. 
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5.  A  Unified  Exposition  of  Some  Classical  Matrix  Theorems 

This  paper  appeared  during  the  period  of  the  current  report:  Linear  and  Multilinear  Algebra, 
(1989),  Vol.  25,  pp.  137-147. 

There  are  several  interesting  and  elegant  theorems  about  hermitian  matrices  that  can  be  unified 
by  a  very  simple  inequality  for  inner  products. 

The  concepts  contained  in  this  paper  are  all  "name"  theorems:  the  Hadamard  determinant 
theorem;  the  Fischer  inequality;  the  Kantorovich  inequality;  Weyl’s  inequalities.  What  may  be  less 
widely  known  is  the  observation  that,  except  for  Weyl’s  inequalities,  these  results  are,  in  fact, 
equivalent  to  one  another,  and  in  turn,  equivalent  to  an  elementary  inequality  usually  referred  to 
as  the  "obvious"  lo..or  bound  in  the  Kantorovich  inequality. 

The  penultimate  section  of  this  note  incorporates  a  brief,  self-contained  discussion  of  compound 
matrices.  This  old  subject  is  less  well  known  than  it  should  be,  and,  as  it  turns  out,  provides 
additional  insight  into  the  theorems  listed  above. 

In  the  final  section,  the  important  inequalities  of  H.  Weyl  that  relate  singular  values  and 
eigenvalues  are  discussed. 

The  statements  of  all  of  the  theorems  appearing  in  the  sequel  are  found  in  [2],  with  details 
concerning  their  origins.  The  analysis  of  equality  in  Weyl’s  inequalities  contained  in  the  last 
section  of  the  paper  are  new. 
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6.  A  Note  on  the  Determinants  and  Eigenvalues  of  Distance  Matrices 

This  paper  appeared  during  the  period  of  the  current  report:  Linear  and  Multilinear  Algebra, 
(1989),  Vol.  25,  pp.  219-230. 

A  set  of  n  vectors  in  a  real  s-dimenaional  Euclidean  space  define  an  n-square  matrix  D  whose  i,j 
entry  is  the  distance  between  the  ith  and  the  jth  vectors.  This  symmetric  matrix  has  0  down  the 
main  diagonal  and  positive  off  diagonal  entries.  Such  matrices  are  called  distance  matrices. 
Distance  matrices  arise  in  a  class  of  techniques  known  as  multidimensional  scaling  (MDS).  The 
purppose  of  MDS  is  to  reconstruct  data  concerning  the  vectors  from  an  examination  of  their 
distance  matrix.  Thus  the  invertibility  of  distance  matrices  (and  related  matrices)  is  an  important 
necessary  condition  on  the  vectors  which  is  directly  discernable  from  their  distance  matrix.  The 
purposes  of  the  present  paper  are:  to  examine  the  behavior  of  the  determinant  of  distance 
matrices:  to  obtain  explicit  formulas  for  the  eigenvalues  of  a  distance  matrix  when  the  vectors  from 
a  regular  polygon  in  real  2-dimensional  space;  and  to  present  in  simplified  from  certain  classical 
results  on  the  rank  and  signatutre  properties  of  matrices  whose  entries  are  squared  distances. 
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7.  Symmetry  Properties  of  Higher  Numerical  Ranges 

This  paper  appeared  during  the  period  of  the  current  report:  Linear  Algebra  and  Appl.,  1988, 
Vol.  104,  pp.  141-164 

Let  A  be  a  linear  operator  on  a  finite  dimensional  unitary  space  V  of  dimension  n.  The  kth  higher 
numerical  range  of  A,  denoted  by  Wk(A),  is  the  totality  of  complex  numbers  tr(PAP)  where  P  runs 
over  all  k-dimensional  orthogonal  projections  on  V.  In  this  paper  it  is  proved  that  Wk(A)  is  a 

polygon  with  the  real  axis  as  a  line  of  symmetry,  k  =  1 . n,  if  and  only  if  A  is  normal  with  a  real 

characteristic  polynomial.  Several  non-normal  examples  are  constructed  in  order  to  investigate 
the  extent  to  which  the  symmetry  of  all  of  the  Wk(A)  is  required. 
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8.  Lower  Bounds  for  the  Norms  of  Decomposable  Symmetrized  Tensors 

This  paper  appeared  during  the  period  of  the  current  report :  Linear  and  Multilinear  Algebra, 
(1989),  Vol.  25,  pp.  269  -  274. 

Lower  bounds  are  given  for  the  difference  of  two  decomposable  symmetrized  tensors.  The  first 
bound  uses  a  norm  which  makes  the  component  vectors  in  a  decomposable  symmetrized  tensor 
part  of  an  orthonomral  basis.  The  second  bound  holds  only  for  decomposable  elements  of 
symmetry  classes  whose  associated  characters  are  linear. 
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